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Adiabaticity and spectral splits in collective neutrino transformations 
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Neutrinos streaming of! a supernova core transform collectively by neutrino-neutrino interactions, 
leading to "spectral splits" where an energy -E sp ij t divides the transformed spectrum sharply into 
parts of almost pure but different flavors. We present a detailed description of the spectral split 
phenomenon which is conceptually and quantitatively understood in an adiabatic treatment of 
neutrino-neutrino effects. Central to this theory is a self-consistency condition in the form of two 
sum rules (integrals over the neutrino spectra that must equal certain conserved quantities). We 
provide explicit analytic and numerical solutions for various neutrino spectra. We introduce the 
concept of the adiabatic reference frame and elaborate on the relative adiabatic evolution. Violating 
adiabaticity leads to the spectral split being "washed out." The sharpness of the split appears to 
be represented by a surprisingly universal function. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

In a dense neutrino gas, neutrino-neutrino refraction 
causes nonlinear flavor oscillation phenomena [1-18]. In 
the region between the neutrino sphere and a radius of 
several hundred kilometers in a core-collapse supernova 
(SN), the flavor content of neutrino fluxes is dramatically 
modified [7-18]. Ordinary MSW resonances typically oc- 
cur at much larger distances and process the spectra fur- 
ther in well-understood ways [19, 20]. 

One particularly intriguing feature of the emerging 
fluxes is a "spectral split." What this means is best 
described with the help of an example. In Fig. 1 we 
show thermal v e and v e flux spectra with an average en- 
ergy of 15 MeV produced in the neutrino sphere of a SN 
core (thin lines). This example is schematic in that we 
assume equal average z/ e and v e energies but an over- 
all v e flux that is only 70% of the v e flux. Moreover, 
the smaller fluxes of the other species v^, v T , Cy, and v T 
are completely ignored. We are considering two-flavor 
oscillations between v e and another flavor v x , driven 
by the atmospheric mass squared difference Am 2 = 2 
3 x 10~ 3 eV 2 and the small 1-3 mixing angle. Therefore 
in Fig. 1 we show the z-components of the usual two- 
flavor polarization vectors, where "up" denotes the elec- 
tron flavor and "down" the x-flavor. In other words, the 
spectrum represents electron (anti)neutrinos where it is 
positive and x (anti)neutrinos where it is negative. 

After a few 100 km, the neutrino-neutrino effects have 
completely died out. The emerging spectra are shown as 
thick lines. The antineutrinos (dashed) have completely 
flipped to the cc-flavor. The same is true for the neu- 
trinos (solid) down to a critical energy -E sp iit where the 
spectrum splits. All v e below this energy emerge in their 
original flavor. Spectral splits of this sort were first ob- 
served in the numerical simulations of Duan et al. [9] 
where a first interpretation was given. Later these au- 



thors have introduced the term "stepwise spectral swap- 
ping" to describe this phenomenon [17]. 

In a previous short paper [14] we have shown that -E sp ut 
is fixed by what we call lepton-number conservation. For 
a sufficiently small in-medium mixing angle the collec- 
tive flavor transformations have the property of conserv- 
ing flavor number in the sense that only pairs v e v e are 
transformed to pairs v x v x , whereas the excess v e flux is 
conserved [11]. While this property does not explain the 
existence of a spectral split, it explains its location that 
for the example of Fig. 1 is at £? sp iit = 9.57 MeV. 

The split phenomenon becomes clearer when the spec- 
tra are represented in terms of the oscillation frequency 
u) = Am 2 /2E instead of the energy itself [14]. In the 
evolution equations the only difference between neutri- 
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FIG. 1: Neutrino spectra at the neutrino sphere (thin lines) 
and beyond the dense-neutrino region (thick lines) for the 
schematic SN model described in the text. Solid: neutri- 
nos. Dashed: antineutrinos. Positive spectrum: electron 
(anti)neutrinos. Negative spectrum: x (anti)neutrinos. 
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FIG. 2: Same as Fig. 1, now with u> = Am 2 /2E as indepen- 
dent variable: co < is for antineutrinos, lo > for neutrinos. 



nos and antineutrinos is that the latter appear with a 
negative lo. Therefore, one should think of the neu- 
trino and antineutrino spectra as a single continuum 
— oo < lo < +oo. The two spectra merge at the point 
lo = that represents both a neutrino or antineutrino of 
infinite energy. In Fig. 2 we show the example of Fig. 1 
in terms of the w-variable where w sp ii t = 0.314 km -1 . 
All modes with lo < w sp nt are transformed whereas the 
ones with lo > w sp nt stay in their original flavor. 

The split occurs among the neutrinos (w sp ii t > 0) be- 
cause a SN core deleptonizes, producing an excess v e flux. 
On the other hand, an excess of the D e -ftux is typical for 
the leptonizing accretion torus of merging neutron stars 
that are probably the central engines of short gamma- 
ray bursts. Here we have w S piit < and the situation is 
reversed: All modes with lo > uj sp ut convert (Sec. Ill C). 
However, we will always use an excess v e flux without 
loss of generality. 

The interpretation of the split phenomenon proposed 
by Duan ct al. [12] motivated us to develop a quan- 
titative theory of this effect based on the idea of adi- 
abatic evolution of individual modes [14]. While the 
self-interacting neutrino ensemble follows a set of non- 
linear equations of motion (EOMs), the evolution of in- 
dividual modes can be written in the form of single- 
particle EOMs where the neutrino-neutrino interaction 
is included self-consistently by a slowly varying "mean 
field" in the single-mode Hamiltonians. The crucial step 
was to obtain two sum rules (integrals over the neutrino 
spectrum in the w-variable) that had to equal certain 
conserved quantities, notably the flavor lepton number 
of the system [14]. These sum rules provide an explicit 
adiabatic solution of the problem. 

The degree of adiabaticity manifests itself in the 
"sharpness" of the spectral split. The adiabatic evolution 
is driven by the decreasing neutrino-neutrino interaction 
strength as a function of distance from the SN core. As 
we dial the scale height to larger values, the spectral split 
becomes sharper, approaching a step function in the limit 
of an infinitely slow evolution. 



This situation is similar to the usual MSW case, for ex- 
ample in the outer layers of a SN envelope. The v e flux 
streaming through a resonance region completely con- 
verts, for the normal hierarchy, while the antineutrinos 
remain unchanged. In other words, in the case of perfect 
adiabaticity the usual MSW effect is a spectral-split phe- 
nomenon with w S piit = 0. In this sense the main effect 
of the neutrino-neutrino interaction is to shift Lo sp ut to a 
non- vanishing value. We will further discuss this analogy 
in the context of adiabaticity breaking in Sec. VII. 

The main purpose of our paper is to elaborate sev- 
eral issues that were kept short in our previous publica- 
tion [14]. In addition, we will pay greater attention to 
the actual definition of what we mean with adiabaticity 
and we will study phenomena caused by deviations from 
a purely adiabatic evolution. We begin in Sec. II with a 
description of the framework for our study: the EOMs 
and our simplifications and approximations. In Sec. Ill 
we present the adiabatic solution in the form of two sum 
rules. The following Sees. IV- VI are dedicated to explicit 
analytic and numerical solutions of the EOMs in partic- 
ular cases. In Sec. VII we study the impact of deviations 
from perfect adiabaticity. We conclude in Sec. VIII. 



II. FRAMEWORK AND SIMPLIFICATIONS 
A. Equations of motion (EOMs) 

Nonlinear flavor transformations caused by neutrino- 
neutrino interactions involve two main complications. 
One is the energy dependence of the effect, a question 
that we will address here. The other is the "multi-angle" 
nature of the transformations. If a neutrino gas is not 
perfectly isotropic, neutrinos moving in different direc- 
tions feel a different refractive effect caused by the other 
neutrinos because the interaction energy between two 
neutrinos depends on their relative direction as (1— cos#). 
As a consequence, one expects kinematical decoherence 
of the flavor content of different angular modes. For 
equal densities of neutrinos and antineutrinos, decoher- 
ence indeed occurs and is even self-accelerating in that 
an infinitesimally small anisotropy is enough to trigger 
an exponential runaway towards complete flavor equilib- 
rium [13]. On the other hand, if the neutrino-neutrino in- 
teraction decreases slowly and if the asymmetry between 
v e and v e is sufficiently large, one observes numerically 
that decoherence is suppressed. In the SN context, the 
deleptonization flux appears to be large enough to pre- 
vent multi-angle decoherence [15]. 

Therefore, we limit ourselves to an isotropic gas with 
a slowly decreasing density. The ensemble is fully char- 
acterized by a matrix of densities in flavor space for each 
spectral component E, separately for neutrinos and an- 
tineutrinos. The energies enter only through the vacuum 
oscillation frequencies lj = Am 2 /2E which determine the 
behavior of different modes. Therefore, it is more conve- 
nient to use lo directly to label the modes. In the two- 
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flavor context we use the usual polarization vectors so 
that the ensemble is represented by two sets of polariza- 
tion vectors P w and P w . The global polarization vectors, 

r-OO r-OO 

P = / dujP^ and P = / do>P w , (1) 
Jo Jo 

are initially normalized as 

P=|P| = 1 and P=|P|=a<l. (2) 

Notice that P w and P w are densities relative to the ui 
variable so that they have the dimension of The in- 
dividual lengths P w = | Pa; | and P w = |P W | are conserved 
as long as oscillations are the only form of evolution, i.e., 
in the absence of dynamical decoherence caused by col- 
lisions. In contrast, the total lengths P and P are only 
conserved when all individual polarization vectors evolve 
in the same way. P^ and P u are the spectra of frequencies 
given by an initial condition. 

We assume that initially all P w and P w are aligned. 
In other words, all neutrinos and antincutrinos begin in 
the electron flavor. It is straightforward to generalize 
this assumption to a situation where some modes begin 
in the x-flavor ("spectral cross-over," see Sec. IV C). In 
this more general case, the P w and P w may initially also 
be anti-aligned relative to the common flavor direction. 

Following the notation of Ref. [11] let us introduce the 
difference vector 

D = P P, (3) 

representing the net lepton number. Using the conven- 
tions of Ref. [21], the EOMs are 

d t Pu> = (+cjB + AL + /iD) x P w , 

dtPu, = (-wB + AL + /iD) x P w , (4) 

where B = (sin 29, 0, cos 29) is the mass direction and 
9 the vacuum mixing angle. For the inverted mass hi- 
erarchy 9 r~j 71-/2. The unit vector L = (0,0,1) is 
the weak charged-current interaction direction so that 
B L = cos 29. Finally, A = \/2G_Fn e is the usual matter 
potential and \x = \[2GFn v the potential due to neutrino- 
neutrino interactions with n e and n v being the electron 
and neutrino densities. We use time t instead of radius 
r because we study an isotropic system that evolves in 
time instead of a single-angle system evolving in space. 

In the SN context, the neutrino flux dilutes as r~ 2 . 
Moreover, the trajectories become more collinear with 
distance, leading to an average suppression of the inter- 
action energy by another factor r -2 [4] so that altogether 
we should use n(t) cx t~ 4 . However, in the adiabatic limit 
the exact form of n(t) is irrelevant as long as its rate of 
change is small. For our numerical examples it is more 
convenient to use the profile 

= no exp(-i/r) . (5) 

The exponential form is simpler because it involves the 
same scale \d\n /i(t) / dt]^ 1 = r at all times. Adjusting 



t allows one to control the adiabaticity of the transi- 
tion with a single parameter. Assuming the atmospheric 
Am 2 and typical SN neutrino energies, we find that 
luq = 0.3 km -1 is typical for the average vacuum os- 
cillation frequency. The crucial action takes place at the 
relatively large radius of a few hundred kilometers. To- 
gether with the r -4 scaling this implies a characteristic 
rate of change of 0.01 km -1 . Therefore, in Eq. (5) a value 
r = (0.03 wo) -1 roughly captures the situation in a SN. 



B. Removing the usual matter effect 

Unless there is an MSW resonance in the dense-neu- 
trino region, one can eliminate AL from Eq. (4) by going 
into a rotating frame [8, 11]. The price for this transfor- 
mation is that in the new frame B has a static component 
P>l = B cos 29 along the L direction and a fast-rotating 
one in the transverse direction. Assuming that the ef- 
fect of the fast-rotating component averages to zero, we 
may consider the equivalent situation where the usual 
matter is completely absent and initially all polarization 
vectors are aligned with B (which really is the projection 
of the true B on the weak-interaction direction L). This 
is equivalent to a vanishing vacuum mixing angle. 

Of course, this is not the complete story because for a 
strictly vanishing vacuum mixing angle, no transforma- 
tion effects would take place at all. The rotating compo- 
nent of B disturbs the perfect alignment of P w and P w , 
providing a kick start for flavor conversions [8]. For a 
system consisting of equal numbers of neutrinos and an- 
tineutrinos, the equations were solved analytically [11]. 
The net effect was found to be equivalent to a small (but 
nonvanishing) effective mixing angle, although the evo- 
lution of the polarization vectors is initially modulated 
by the rotating B field with the frequency A. We have 
here a different system where the density of neutrinos ex- 
ceeds that of antineutrinos and where the evolution of the 
system is driven by the fi(t) variation. Numerically one 
finds that once more the effect of matter is equivalent to 
reducing the vacuum mixing angle to a certain effective 
mixing angle 6 c ff similar to the usual in-medium mixing 
angle [15], but an analytic treatment of how ordinary 
matter modifies the initial evolution is not available. 

Motivated by this evidence we make the simplification 
of ignoring the ordinary matter term entirely, schemati- 
cally accounting for it by a small initial effective mixing 
angle. Notice that in the absence of the interaction vec- 
tor L we define the effective mixing angle as the angle 
between D and B. Our treatment is exact for a system 
where the only interaction is that of the neutrinos among 
each other. While we cannot claim with mathematical 
rigor that this system is a faithful proxy for what hap- 
pens in a real SN, numerical observations suggest that 
this is the case. 
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C. EOMs in the w variable 

The equation of motion for P w differs from that for 
P w only by the sign of frequency: u> — > — w. So, the 
only difference for antineutrinos is that in vacuum they 
oscillate "the other way round." Therefore, instead of 
using P w we may extend P w to negative frequencies such 
that P w = P -u (oj > 0) and use only P w with — oo < 
uj < +00. After eliminating the matter term AL the 
EOMs then take on the simple form 

P^ = (wB + /iD) x P w . (6) 

A compact expression for the difference vector is 



/ + OC 
-00 



dusu P w , 



(7) 



where 



s w = sign(w) 



10 



(8) 



Integrating both sides of Eq. (6) over s^duj provides 



where 



D = B x M. 



f +00 

M = / duj s^uj P u 



-f 

J — ( 



(9) 



(10) 



is the "effective magnetic moment" of the system. 

The EOM for D is not a closed differential equation. 
However, when fj, is large, all P w (—00 < u> < 00) remain 
pinned to each other, and therefore M oc D. In this case, 
according to Eq. (9) the collective vector D precesses 
around B with the synchronization frequency [6] 



_ M _ 

^synch — — 



f 00 dto i 
j —00 



(11) 



If /i is suddenly turned off, all P w henceforth precess 
with their individual uj. The transverse component of 
D quickly averages to zero (kinematical decoherence) , 
whereas the component along B is conserved. We will 
here study the opposite limit where \i decreases slowly, 
leading to a very different final result. 



D. Lepton-number conservation 

Since B is constant, Eq. (9) shows that d t (D ■ B) =0 
and therefore 



D\\ = B D = const., 



(12) 



where here and henceforth we will denote the component 
of a vector parallel to B with the index ||, the transverse 
component with the index _L. 



The physical interpretation of the conservation law 
Eq. (12) is what we call "flavor lepton-number conser- 
vation" or simply "lepton-number conservation." Rigor- 
ously it only means that in the absence of matter the 
projection of D on B is conserved. However, in all exam- 
ples of practical interest in the SN context, the effective 
mixing angle is small so that initially D and B are nearly 
collincar. Therefore, the net v e flux from delcptonization 
is approximately conserved until an MSW resonance is 
encountered. Collective effects only induce pair transfor- 



mations of the form v e v e 



v x v T , whereas the excess of 



the v e flux from deleptonization is conserved. 



III. ADIABATIC SOLUTION 

A. Adiabaticity and co-rotating plane 

To understand the evolution of our system in the limit 
of a slowly changing fi, we rewrite the EOMs as 



(13) 



where —00 < u < +00. We have introduced an "individ- 
ual Hamiltonian" for each mode 



H w = wB + /iD . 



(14) 



According to the standard notion, the evolution is adi- 
abatic when the rate of change w off (H w ) of each H w is 
slow compared to the precession speed w cfr (P w ) of P w : 



w eff (P w )»u; ctt (H w ). 



off/ 



(15) 



In this case P w follows H w . In general P w need not 
coincide with H w but rather moves around H w on the 
surface of a cone whose axis coincides with H w and the 
cone angle is constant. 

Let us check the adiabaticity condition for our system. 
Since \D\ < 1 + a, the speed (effective frequency) of the 
P w motion can be estimated according to Eq. (13) as 



off 



(P„) < ^ 2 +^ 2 (l + a) 2 - 



(16) 



D and therefore H w process around B with the synchro- 
nization frequency of Eq. (11): 



UJ cii (H w ) — ^synch ~ ^Oj 



(17) 



where u>o is some typical frequency in the spectrum. The 
Hamiltonian H w changes both because of the /^-decrease 
and because of the precession of D. For \i <~ uj or smaller 
we have 



^ ff (P w )^c ott (H w ), 



off/ 



(18) 



and therefore the adiabaticy condition that requires 
much faster motion of P w is not fulfilled. Apparently 
though the adiabaticity condition is fulfilled for /i 3> loq- 
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However, it is usually not stressed that the question of 
adiabaticity depends on the coordinate system and there- 
fore adiabaticity has a relative meaning. For a single po- 
larization vector the evolution is trivially adiabatic in a 
reference frame co-moving with its evolution, although 
this observation is useless as a means for solving the 
EOM. It is also trivial to find a frame where the evo- 
lution is non-adiabatic. Even for the case of the ordinary 
MSW effect we can go into a frame that rotates with 
a large frequency around the z-direction. In this frame 
the evolution does not look adiabatic. Of course, the fi- 
nal outcome is unchanged because the physical evolution 
does not depend on the coordinate system in which it is 
discussed, but the interpretation is different. In the fast- 
rotating frame the ordinary MSW effect is interpreted 
as a parametric resonance caused by a B-field with a 
fast-rotating transverse component. 

Let us return to the case at hand. Inspecting Eq. (14) 
reveals that the fast motion of the single-mode Hamil- 
tonians in the "laboratory frame" is caused by the fast 
motion of the single common vector D. If the motion of 
D is essentially a precession around B as suggested by 
Eq. (9), the relevant co- moving frame is simply the plane 
of B and D that rotates around B. In other words, while 
both H w and P w precess fast around B, in the co-rotating 
frame they may move slowly relative to each other and 
relative to the common frame which then plays the role 
of the adiabatic frame for all P w . 

Going to a rotating frame to simplify the EOMs was 
first suggested by Duan et al. [8], but later they have 
questioned that the evolution can then be called adi- 
abatic because it is not adiabatic in the laboratory 
frame [16]. We believe that the notion of adiabaticity is 
nevertheless appropriate because the laboratory frame is 
not special. (In fact, we have already transformed away 
the usual matter effect by going into a rotating frame 
which for us plays the role of the laboratory frame.) The 
evolution of the multi-mode system should be called adi- 
abatic if there exists a common frame in which all modes 
evolve adiabatically. 

In the present context, the co- rotating frame spanned 
by B and D naturally provides a common adiabatic 
frame. Independently of the motion of the individual P w 
that define the overall D, the single- mode Hamiltonians 
H w always lie in a single plane spanned by the vectors B 
and D. Relative to the laboratory frame, this common 
or co-rotating plane moves around B with the instanta- 
neous co- rotating frequency u c . The Hamiltonians in the 
co-rotating frame are 

H w = (w - lo c ) B + fiD . (19) 

We use the same notation in both frames because the 
relevant components H u », H u ±, D», and D± remain in- 
variant under transition from the one frame to another. 

Let us find the co- rotating frequency uo c . By assump- 
tion all P^, and consequently M, stay in the co-rotating 
plane. Therefore we can decompose 



and rewrite the EOM of Eq. (9) as 
d t T> = lu c B x D. 



(21) 



Projecting Eq. (20) on the transverse plane we find uj c 
M±/D± or explicitly 



r°°du 

J — oo 



/-TO diJ s " p ^ 



f + °° cLos w cjP„a 



(22) 



When fi — > oo and all P w are aligned, this is identical 
with the synchronization frequency Eq. (11). 



B. Adiabatic solution in terms of sum rules 

Initially when ji is very large, all individual Hamil- 
tonians are essentially aligned with D. In turn, D is 
aligned with the weak-interaction direction if initially all 
polarization vectors P w are aligned with that direction. 
In other words, all neutrinos are prepared in interaction 
eigenstates and initially P w oc H w . The adiabatic evo- 
lution would imply that if \i changes slowly enough, P^ 
follows H w (fj.) and therefore remains aligned with H w (fi) 
at later times as well. So the adiabatic solution of the 
EOMs for our initial condition is given by 



P»M = H u (/i)i>, 



(23) 



where P w = |P W | and = H W /|H W | is a unit vector in 
the direction of the Hamiltonian. 

According to the solution Eq. (23) all P w being con- 
fined to the co-rotating plane evolve in this frame accord- 
ing to the change of /U. The solution is implicit since H w 
depends via D = D(P W ) on the polarization vectors: 



H W (P W ), 



P W (H W ). 



(24) 



We will see that there is a consistent solution which sat- 
isfies these relations. Notice that by imposing the "align- 
ment" condition Eq. (23) on P w , we also restrict the 
Hamiltonian, since H w depends on P w . 

Let us find explicit solutions P u « and P u ±. Since P w is 
conserved and simply the spectrum of neutrinos given by 
the initial condition, we only need to find the Hamiltoni- 
ans Eq. (19). They, in turn, are completely determined 
by uj c (fi) and D±(/j,) because D^y is known to be constant 
as determined by the lepton number. In other words, we 
need to find self-consistent results for w c {v) an d D±(^i), 
the component transverse to B, for a given Du. 

The implicit adiabatic solution Eq. (23) gives us two 
equations for each u> which we use to find u) c (/j,) and 
D±(fi). Projecting Eq. (23) onto the perpendicular and 
parallel directions with respect to B provides 



M = bB + w c D 



(20) 



Pu>± 



p»\\ = 



H. 



(25) 
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From Eq. (19) wc infer 



so that 



^11 = 



(u -u> c + fiD ll )P LO 
^/(uj-uj. + fiD^ + ipD^ 



y/(w-u c + nD\\)* + (jiD ± )* ' 
Integration of these equations over s^duj gives us 

{w-w c + iiD\\)P u 



+ 00 



duj s u 

i 

duj s u 



P 

r UJ 



(26) 

(27) 
(28) 

. (29) 



^/[(u,-co c )/^ + D^+Dl 



• (30) 



For a given Du and spectrum P w , we can determine uj c 
and D± from Eqs. (29) and (30) for any \x and thus find 
explicit solutions. 

Equations (29) and (30) can be considered as sum 
rules — integrals over frequencies that should be equal to 
certain conserved numbers. They can be added to each 
other with various factors so that one can rewrite them in 
different forms depending on convenience. For example, 
Eq. (29) can be rewritten with the use of Eq. (30) as 



r + oo 

Ld r = I duj S u 



-J 

J —I 



LUP„ 



[UJ — Ul, 



(31) 



which coincides with the form derived in our previous 
paper [14]. 

While different forms of the sum rules are equivalent, 
their limiting behavior for fj, — > or /! — > oo is not always 
equally manifest. For example, the physically intuitive 
expression Eq. (22) for uj c is valid only for D±(fi) ^ 
while for D± = it has a 0/0 feature. On the other hand, 
Eq. (31) for fi — > oo is equivalent to the synchronization 
frequency Eq. (11) without any ambiguity. 

Inserting the solution Eq. (23) into the EOM Eq. (13) 
we obtain 



d t ujp) = o 



(32) 



which is satisfied exactly if d t ^ = and <9 t D = 0. That is, 
the medium has constant neutrino density and the differ- 
ence vector does not change. The latter is satisfied when 
all polarization vectors stay unchanged, but this may not 
be the only possibility. We will call this self-consistent so- 
lution the "static solution." For a slowly varying density 
the condition Eq. (32) is only approximately satisfied. 
As in the usual MSW case, Eq. (23) provides an approx- 
imate solution of the EOMs, when the dependence of the 
Hamiltonian on time is negligible. The adiabatic solution 



(the lowest order term in the adiabatic perturbation the- 
ory) is given by a continuous set of static (instantaneous) 
solutions. 

Let us comment on the static solution. It consists of 
all P w being confined to the co-rotating plane. The mo- 
tion of all polarization vectors is limited to a common 
precession around B with frequency ui c , each of them 
having a different zenith angle relative to B. Duan et al. 
have termed this form of motion the "pure precession 
mode" [12]. In contrast to the usual MSW case, the 
static solution implies an additional condition D = const, 
which does not mean in general that the individual polar- 
ization vectors are static. Another issue is the stability 
of the solution since D is a dynamical system and not 
a rigid vector. Furthermore, separate equalities dtn = 
and d t T> — may not be the only solutions of Eq. (32). 
Our numerical studies show, however, that the static so- 
lutions are not only self-consistent, but also stable. 

For every spectrum P w the pure precession solutions 
are characterized by the two parameters D\\ and fi. For 
fj, — > oo all Hamiltonians are aligned and thus all P w 
are aligned as well, provided that they are in a pure pre- 
cession mode. Assuming the neutrinos are produced in 
interaction eigenstates, all P w are tilted relative to B 
with twice the effective mixing angle so that 



r+oo 

D H =coa29 cf i I dus^Poj 



(33) 



In this sense we have, for every spectrum P u , a two- 
parameter family of static solutions given in terms of 
cos2# ff and \i or cquivalcntly, in terms of Dn and /i. 



C. Split frequency 

In the limit [i — > oo all polarization vectors are aligned 
with each other in a direction given by the choice of 
cos2#^j., i.e., by the initial condition. In the opposite 
limit, jjb — * 0, the solution is given by 



IT 



(34) 



where wj? = w c (m ~* 0). All Hamiltonians and thus all 
P w with lo > lo c are aligned with B, whereas those with 
uj < w° are anti-aligned. Therefore, we have a spectral 
split at the frequency 



W sp lit 



(35) 



which usually is not equal to zero. 

D± = because in the end all polarization vectors 
are (anti)aligned with B. Moreover, for fi — > Eq. (29) 
approaches a well-defined limit, 



Dn 



= / 
= / 



(U) - W S plit) 

duj s w p = Pu, 

^{uj- w S piitr 



(36) 
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For our usual case of an excess neutrino flux (D\\ > and 
w sp iit > 0) this equation is identical with 

/0 (•"split r + ao 

Pudw- P L ,du + P u du . (37) 

-OO JO •'"split 

The first term enters with a positive sign because antincu- 
trinos enter lepton number with a negative sign, but are 
now oriented opposite to B, the next term has a nega- 
tive sign because it represents those neutrino modes that 
in the end are oriented against B, and the third term 
has a positive sign because it represents neutrinos with 
orientation along B. 

In practice we are mostly interested in the application 
to SN physics where the effective mixing angle is small. 
In this case there are other intuitive ways to state the 
condition of lepton-number conservation and to deter- 
mine the split frequency. We can rewrite the definition 
of D in Eq. (7) for 9^ w in the following explicit form 



Instead of Eq. (40) we get 



/0 /•"split r 

Pu>duj+ P ul dcu + 

-oo J J oj 



+ oo 

1\ - / PujdLU. (38) 

"split 



Now we have two equations, Eq. (37) and Eq. (38), for 
integrals over different parts of spectra. The first equa- 
tion gives us the total lepton number of the initial state 
where antineutrinos enter with a negative sign, neutrinos 
with a positive sign. The second equation simply reflects 
that all modes below o-> sp iit were flipped. 

Subtracting or adding these equations eliminates dif- 
ferent parts of the spectrum and reveals 



/0 /-"split 
dLoP u = 
-oo J0 



duj P u (39) 



and 



f + 00 

Du = I dcuP^ 



(40) 



Equation (39) has the interpretation that all antineutri- 
nos (— oo < u < 0) are converted and therefore to con- 
serve lepton number the same number of neutrinos should 
be converted. From the continuity in u> we infer that 
those are neutrinos with frequencies with < u < u> sp nt- 
Equation (40) has the interpretation that the net lepton 
charge of the system should be determined by the high 
frequency part of the neutrino spectrum. 

For the case of a dominant antineutrino flux (a > 1) 
we have D» < 0, the split frequency is negative and the 
spectral split occurs in the antineutrino channel. Now 
all modes with frequencies above w sp i;t are converted 
whereas the modes with negative u below Lo sp \\ t stay un- 
changed. The sum rule for the determination of w sp iit 
now becomes 

/"split r0 r+oo 

dcuP LJ + dwP u - dwP^. (41) 
-OO •'"split J^O 



f 



"split 



divPu, . 



(42) 



The net charge now comes from negative frequencies. So 
in this case the entire neutrino flux and high energy part 
of the antineutrino flux change flavor whereas the low 
energy antineutrino flux is unchanged. 



IV. SOLUTIONS FOR NEUTRINOS 
A. Box spectrum 

As a first example for explicit adiabatic solutions 
we consider the simplest possible case where neutrino- 
neutrino interactions lead to a spectral split in the adi- 
abatic limit: An ensemble consisting of neutrinos only 
(a = 0), no antineutrinos and no ordinary matter. Now 
D = P and correspondingly D\\ = P|| and D± = P±. 
Of course, this is an abstract example meant to illustrate 
the general phenomenon. In subsequent sections we will 
turn to realistic situations relevant in a SN. 

We use a flat distribution of oscillation frequencies with 
the average u>o, i.e., the box-like spectrum 



P" 



f (2w )" 
1 1 







for < oj < 2uj , 
otherwise. 



(43) 



With Py the conserved projection of P on the B-direction 
we find from Eqs. (31) and (37), recalling that Eqs. (39) 
and (40) are valid for 6£g w 0, 



ojo x 



1 

(1 



for [i 
for \i 



> oo, 
0. 



(44) 



Notice that for P\\ = 0, corresponding to maximal initial 
mixing, the frequency ui c = ojo is constant for the en- 
tire evolution. For Pn = 1 (vanishing mixing angle) we 
have lu c — 0, corresponding to the absence of any flavor 
evolution. We will see in Sec. V that in the presence of 
antineutrinos even a very small mixing angle leads to a 
large effect, whereas here we need to assume a non-zero 
mixing angle to obtain something visible. 

In Fig. 3 we show the initial and final distribution P^y 
for the example Py = 0.5, so that according to Eq. (44) 
^spiit = = 0.5a;o. The dotted line indicates the fi- 
nal distribution in the adiabatic limit, the solid line is 
the result of a numerical example with /j,(t) from Eq. (5) 
with r = (0.03 wo) -1 - The spectral split indeed becomes 
sharper with increasing r (increasing adiabaticity) as dis- 
cussed in more detail in Sec. VII. 

We may also follow the evolution of the individual po- 
larization vectors. In Fig. 4 we show P w || for 51 equally 
spaced modes as a function of uj /h for the example of 
Fig. 3. At first, all of them begin with the common start- 
ing value P w || = 0.5. Later they spread out and eventu- 
ally split, some of them approaching +1 and the others 



8 



1 - 



c 

o 



o 

* 0.5 



Q. 
I 

CD 



- 



3 -0.5 - 



o 

Q. _1 



-i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 



J I I I I I I I I I I I I I I I I L 



0.5 1 1.5 2 

FIG. 3: Initial (thin line) and final (thick line) neutrino spec- 
tra for the example of Eq. (43). Dotted: fully adiabatic. Solid: 
numerical example with r = (0.03wo) _1 in Eq. (5). 
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FIG. 4: Evolution of 51 equally spaced modes for the same 
example as in Fig. 3. Top: Numerical solution with r = 
(0.03 wo) - • Bottom: Adiabatic solution derived in Sec. IV B. 



— 1 as predicted. Some of them first move down and then 
turn around as a result of u c changing as a function of 
time. As a result of imperfect adiabaticity, a few modes 
do not reach ±1 but get frozen before this destination. 



B. Analytic solution 

In the adiabatic limit, the same results follow from our 
explicit solution of the EOMs. For the box spectrum 
Eq. (43) the integrals Eqs. (29) and (30) are easily per- 



formed. After some transformations we find 

P± 



e K — e~ 
2k 



where 



K = iVo/fl. 



(45) 



(46) 



For fi — > oo and fi — > the limits of lo c agree with the 
results predicted in Eq. (44). For /j^oowe have P± = 
(1 — P 2 ) 1 / 2 , representing the initial condition P = 1. For 
/i^Owc find P± = 0, corresponding to all polarization 
vectors either aligned or anti-aligned with B. 

For our example with P|| = 0.5 we show lu c (k) and 
P_i_(k) in Fig. 5. According to Eq. (45), iv c decreases 
monotonically with increasing k from uj c = loq to lu c — 
w S piit = 1 — P\\ = 0.5. The transverse component P± 
decreases from (1 - Pf) 1/2 = V075 = 0.866 down to 0. 
With Eq. (27) the analytic results provide exact adiabatic 
solutions for the || -component of the polarization- vector 
spectrum. In the bottom panel of Fig. 4 we show the 
evolution for 51 equally spaced modes in analogy to the 
numerical result shown in the upper panel. The agree- 
ment for not very small /x is striking, thus confirming the 
correctness of our picture based on the adiabatic evolu- 
tion in the co-rotating frame. The curves obtained in the 
adiabatic approximation and by numerical solution of the 
original evolution equations do not agree for the modes 



1 — I I I Mill 1 — I I I I I III 1 — I I I III 
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FIG. 5: Analytic results for ui c (n) and P±(k) according to 
Eq. (45) with the initial condition Pn =0.5. 



9 



close to the split (to w w sp nt) at low neutrino densities 
(/i < Wo) where the evolution becomes non-adiabatic. 

Figure 4 shows that the trajectories of different modes 
literally split in the course of evolution, giving us another 
(dynamical) motivation for the term "spectral split" to 
describe the phenomenon. 

It is remarkable that the split is essentially complete 
only for ujq/h 7 or [i < 0.15 wo- In other words, the 
overall motion remains collective in the sense of all polar- 
ization vectors precessing in a common plane for a sur- 
prisingly small value of /z, even for imperfect adiabaticity 
as in the upper panel of Fig. 4. 

Actually for any value of jj,, no matter how small, our 
adiabatic case is an exact solution of the EOMs where 
all polarization vectors precess in a common plane. Of 
course, for very small \x their angles relative to the posi- 
tive or negative B direction are exponentially small, but 
conceptually this is still a collective motion. In the per- 
fectly adiabatic case, collectivity is never lost. Instead, 
the polarization vectors orient themselves in the positive 
or negative B-direction. 

If \x were to switch off suddenly, the transverse com- 
ponent of the overall P disappears by kinematical deco- 
herence because all P w precess with different frequencies. 
In the perfectly adiabatic case, the transverse component 
disappears because all individual P w (anti) align them- 
selves with B while remaining in a common plane. There 
is never any kinematical decoherence caused by differ- 
ences of precession frequencies between different modes. 
The persistence of collectivity for arbitrarily small [i is 
the most surprising aspect of this system. 



C. Spectral cross-over 

At the neutrino sphere in a SN it is expected that the v e 
flux dominates for low energies and the v x flux for high 
energies. This would be the case, for example, if both 
species are emitted with equal luminosities but different 
average energies. In this case the initial P w are cither 
aligned or anti-aligned with the flavor direction. In the 
previous discussion we must then be careful about abso- 
lute signs of the various projections of the polarization 
vectors. 

An important difference occurs in Eq. (23) because 
P w can be anti-aligned with H w so that P w should be 
taken negative for those modes that initially contain the 
"wrong" flavor. The final state of these modes will be 
opposite from where they would have ended if they had 
been occupied with the "correct" flavor. 

The final spectrum is again determined by flavor 
lepton-number conservation. Imagine that the low-w 
modes are initially in the other flavor, corresponding to a 
fraction £ of all modes, and assume a box spectrum (thin 
line in the top left panel of Fig. 6). These modes must 
now end up at P u \\ — +1, affecting the net flavor lepton 
number. The modes above Co> S piit also must end up at +1, 
whereas intermediate ones end up at —1. It is now easy 




« /m « /m 

FIG. 6: Examples for crossed-over neutrino spectra. Top: 
Thin solid: initial. Thick dotted: adiabatic final. Thick solid: 
numerical final with r = (0.03o;o) _1 in Eq. (5). Bottom: 
Numerical evolution of 51 equally spaced modes. 

to work out that in the end (ji — > 0) 

[1 - cos/3 + 2£(1 + cos /?)], (47) 

where (3 is the initial angle between B and P. For £ = 0.2 
and cos/3 = 0.5 we thus find w sp iit = 1.1 wo- In this way 
we predict for the example of Fig. 6 (left column) the final 
spectrum indicated by a thick dotted line. Numerically 
we find the thick solid line. 

The split at w sp iit is somewhat washed out because the 
evolution is not fully adiabatic. On the other hand, the 
transition at the initial spectral cross over remains sharp 
because the end state of those modes depends only on the 
initial sign of P u , but otherwise we are in an adiabatic 
regime of the spectrum. 

In the bottom panel we show the familiar pattern of the 
P u « evolution with the novel feature that the modes that 
were initially in the other flavor cross over to P^ = +1. 

In the right-hand panels of Fig. 6 we show an analogous 
example where the high-u> part of the initial spectrum is 
crossed over. The explanation is analogous. 



V. NEUTRINOS AND ANTINEUTRINOS 
A. Solution for the double-box spectrum 

As a next case towards a more complete understand- 
ing we include antineutrinos in our model, upholding the 
condition that all polarization vectors P w , now with fre- 
quencies in the range — oo < uj < +oo, initially point 
in the same direction. We assume that the fraction of 
antineutrinos is a < 1 of the neutrinos. We will always 
assume the inverted mass hierarchy so that a large flavor 
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transformation effect is caused by the instability of the 
inverted flavor pendulum [11, 16]. 

As a generic example we extend the box spectrum of 
the previous section to include antineutrinos. To be spe- 
cific, we assume a double-box spectrum of equal width, 
but different height for the neutrinos and antineutrinos: 

I C a for — 2(jj < w < 0, 
Pu = y~ x \ 1 for < w < 2uj , (48) 
1 otherwise. 

We will assume a — 0.7 as a standard case. 

The thin line in Fig. 7 represents this initial spectrum 
where antineutrinos are represented with negative fre- 
quencies. Using the same arguments as in the neutrino- 
only case, in the co-rotating frame the antineutrinos are 
modes with even more negative frequencies than those 
neutrinos with u < w sp iit- Therefore, they should simply 
reverse their direction. This is borne out by the numeri- 
cal final spectra in Fig. 7 (thick solid line). The adiabatic 
limit (thick dotted line) is again explained by the conser- 
vation of flavor lepton number. 

Next we consider the evolution of the individual modes. 
In the upper panels of Fig. 8 we show P u h (ii) for neutrinos 
and in the lower panel for antineutrinos. As expected, 
they evolve at first collectively and both global polar- 
ization vectors tilt against the B direction, the familiar 
result of neutrino-neutrino refraction in the inverted hi- 
erarchy case. Later the modes separate to produce the 
final spectrum. In Fig. 8 we juxtapose the direct nu- 
merical solution of the EOMs, assuming sin2i9£§- = 0.05 
and a not very adiabatic r = (0.1 ujq)^ 1 , with the fully 
adiabatic solution for sin 2$^ = (right panels). 

We could have illustrated this case by using the same 
length for the P w for both neutrinos and antineutrinos, 
but occupying only 70% of the antincutrino modes. In 
this alternative double-box case we could have shown the 
evolution of all modes in a single plot because the an- 
tineutrino modes seamlessly join the neutrino modes. 

A new feature of the left panels of Fig. 8 is the "wig- 




-2-10 1 2 

u/u 

FIG. 7: Neutrino spectra for an initial box spectrum with 
70% antineutrinos and sin 29^ = 0.05. Negative frequencies 
correspond to antineutrinos. Thin line: initial. Thick dotted: 
final adiabatic. Thick solid: numerical example with r = 
(0.1 coo)- 1 in Eq. (5). 
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FIG. 8: P w ||(/x) for individual modes for the case of neutrinos 
plus antineutrinos. Left: Numerical solution from the EOMs 
as in our previous paper [14] for sin2(9^f — 0.05. Right: Adi- 
abatic solution for sin2#^ = 0. In each case neutrinos with 
51 modes (top) and antineutrinos with 6 modes (bottom). 



gles" in the curves that stem from the nutation of the 
global flavor pendulum. We will return to this issue later 
but for now only note that the wiggles disappear in the 
numerical solution if we take the scale parameter r very 
large, much larger than is realistic in the SN context. In 
this case the split would be extremely sharp. In order 
to show a not fully adiabatic solution we have here used 
parameters that do not have too many nutation periods 
(relatively small r) and where the nutation depth is not 
too large (mixing angle not too small) or else the plot 
would be too cluttered. In any case, there are no nuta- 
tions in an "extremely adiabatic" situation which for the 
moment we focus on. 

The adiabatic solution for the individual modes shown 
in the right-hand panels of Fig. 8 were found from solving 
our sum rules Eqs. (29) and (30) numerically. The curves 
w c (/i) and D±(fj.) for a vanishing mixing angle are shown 
as thick solid lines in Fig. 9. It shows interesting new 
features compared to the neutrino-only case. For /i — > oo 
we have D± = due to the assumption of a vanishing 
mixing angle. The precession frequency corresponds to 
the expected synchronization frequency However, for 
decreasing \i the precession frequency increases up to an 
apparent cusp, and then decreases eventually down to 
the expected split frequency. At the cusp, D± begins to 
deviate from zero whereas D± = for < < aC 1 ^. 
In other words, the cusp marks the neutrino-neutrino 
interaction strength where the polarization vectors begin 
to tilt away from the B direction and begin to spread in 
the zenith-angle direction. 

For comparison we show, as a thin line, the same result 
for a large initial mixing angle, sin 20£| = 0.5, implying 
that initially D± = (1 - a) sin 20£§ = 0.15. The curve 
is qualitatively similar but does not show any cusps or 
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kinks. Notice that for lo c > 2ujq the co-rotation frequency 
is above the maximal frequency of the spectrum for k < 
1. Therefore, a true split of the modes (in contrast to a 
zenith- angle spread) only starts when uj c = 2uj , that is 
at k = 1. Apparently the dependences in Fig. 9 differ 
from that in Fig. 5 for pure neutrino case. 



Practically it is not possible to extract D± and lu c as 
functions of /i in closed form for generic values of param- 
eters. 

In a SN, the most interesting case is sin 29^ w 0, where 
Dh = 1 — a. This extreme case entails a number of 
simplifications, but also a number of subtleties. It is easy 
to find uj c (k) for special values of k. 
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FIG. 9: Numerical solution of Eqs. (29) and (30) for the 
double-box spectrum of neutrinos plus antineutrinos with 
a = 0.7 and sin 2#2| = 0. The vertical dotted line in- 
dicates Kcusp = 0.013340 below which the system is per- 
fectly synchronized. The horizontal dotted lines indicate 
^synch = = (17/3) wo and w S piit = w° = (14/10) w . The 
thin solid line is the same case with sin 20£§ = 0.5. 



1. Synchronization frequency 

From Eq. (11) we easily derive the initial common os- 
cillation frequency, corresponding to the well-known syn- 
chronization frequency, 



Wsynch — ^0 



1 + a 

l-a 



(51) 



For a = this is o%nch = ^ <^o = 5.667ojo > 2ojo- 
Therefore, initially uj c is "outside of the spectral box," 
i.e., all neutrino and antineutrino modes have negative 
frequencies in the co-rotating frame. 



2. Split frequency 

Later the modes split and the final split frequency at 
/i = is found from v e number conservation to be 



Wspiit = w (1 - D z + a) = lu 2a , 



(52) 



where we have used that for a very small mixing angle 
D|| = 1 — a. With a = 0.7 this is o> sp iit = 1.4 wo m 
agreement with Fig. 7. Independently of the choice of 
a, the final w S piit always lies "within the box" because 
ujo 2a < 2uj so that a split always occurs. In the end 
D± = due to the (anti)alignmcnt of all polarization 
vectors with B. 



B. Properties of the pure precession mode 

Before turning to more physical interpretations of 
these features, notably the cusp, we first investigate 
some mathematical properties of these solutions. For the 
double-box spectrum, the sum rules can be evaluated an- 
alytically. Eq. (30) gives 



(» + + ^yl + Dl) („_ + ^yl+D\) 



l+a 



e 2K , (49) 



where y± = yo ± 2k and y = D\\ — (uj c /uj ) k. The 
integration of Eq. (29) leads to 



\jv\ + D\ + a^yl + D\ - (1 + a)y/y^ + D 2 ± = 2kD {1 . 

(50) 



3. Cusp 

A final special quantity is the strength of the neutrino- 
neutrino interaction where the system exits the synchro- 
nization regime and enters the "bipolar regime." This 
transition is infinitely abrupt when sin 2^ = and cor- 
responds to the cusp in u) c (k) shown in Fig. 9. 

For < k < Kcusp we note that D± = solves Eq. (50) 
identically. Equation (49) can be written in the form 

(l-u) a (l + u) =e 2K , (53) 



where 



2k 



Dh - KLO C /U) 



(54) 



At k = K CUS p the function lu c (k) appears to have a ver- 
tical tangent. If this is true, u(k) also should have a 
vertical tangent or k(u) should have a horizontal one. 



12 



Equation (53) is trivially inverted to provide an explicit 
expression for k(u). Using the condition dn/du — we 
find M CUS p- Sticking this back into n{u) we obtain 

2 \ . / 2a 



cusp 



log 



+ a log 



l+aj °\l+a 
0.01330 (for a = 0.7). 



(55) 



The co-rotation frequency at the cusp is found to be 



1 



1 



LO 



^cusp 1 ^ 

= 11.2148 (for a = 0.7). 



(56) 



The values of K cusp and Lu CUS p are in perfect agreement 
with the numerical result shown in Fig. 8, thus confirming 
that in the cusp ui c has a vertical tangent. 

We observe that Eqs. (53) and (54) allow us to find the 
curve uj c (k) analytically in parametric form, i.e., we can 
extract a pair of functions uj c (u) and k{u). One easily 
shows that the function n{uo c ) does not end at the cusp, 
but rather is well defined for w sync h < lj c < oo. Turning 
this around, the function lu c (k) is double valued for < 
k < K cusp . We show the full solution in the upper panel 
of Fig. 10 where the unphysical part is continued as a 
dashed line. 

In addition, we find numerically that the high-K branch 
also has a smooth continuation below the cusp. This 



K 



100 



10 



1 

0.4 



0.3 E- 
0.2 

0.1 E 





cusp 



- 1 V 1 Mill 1 1 

: \\ ■ 


1 1 1 llllj 


1 1 1 llllj 


1 1 MMU 












1 : i 


1 


1 






CO- 



CO 



synch 



CO 



split 




0.001 0.01 0.1 1 10 

K = O /fJ, 

FIG. 10: Numerical solution of the sum rules Eqs. (29) 
and (30) as in Fig. 9, now including the unphysical solutions 
as dashed lines. 



branch is also unphysical, corresponding to solutions with 
D\ < shown as a dashed line in the lower panel of 
Fig. 10. The cusp actually corresponds to a cross-over of 
two solutions, although the physical solution indeed has 
a cusp. 

The unphysical solutions persist for nonvanishing mix- 
ing angles. We show the complete numerical solutions 
for sin20^j = 0.5 as thin lines in Fig. 10. The physical 
{D\ > 0) and unphysical (D 2 ± < 0) solutions are now 
fully separated: there are no kinks or cusps as far as one 
can tell numerically. 

For sin 28^ — the physical solution has a cusp, D± = 
is exact for the entire range < n < K C usp, yet uj c varies 
in this range as a function of n. These surprising features 
are more easily understood in a simpler case where one 
uses only two polarization vectors, one for neutrinos and 
one for antineutrinos, i.e., for a monochromatic neutrino 
spectrum. We will consider this example in Sec. VI and 
return there to a discussion of the u> c behavior. 



VI. TWO POLARIZATION VECTORS 

As a last explicit example we study a system con- 
sisting of only two polarization vectors Pi for neutri- 
nos and P2 for antineutrinos with Pi = \Pi\ = 1 and 
Pi = |F*2 1 = a < 1- This example allows us to make 
a connection to the previous literature. Moreover, this 
system is a useful toy example to understand some of the 
peculiarities of the general neutrino-antineutrino system, 
notably the cusp feature in the common precession fre- 
quency w c (m)- Finally, realistic spectra do not extend 
to infinite energies and therefore do not have a continu- 
ity in the region around w = 0. In this respect the two 
line-spectrum represent reality to a certain extend. 



A. Pure precession mode 

According to Eq. (4) the equations of motion can be 
written in the form 



Pi = +UJ Q B X P x + /iPi X P 2 , 
P 2 = -LU B X P 2 + /iPi X P 2 . 



(57) 



If Pi and P 2 were free, both of them would precess with 
frequency coo, but in opposite directions. 

We first investigate the "pure precession mode" [12] 
of this system, at first without reference to adiabaticity, 
where both polarization vectors precess in a single plane 
around B. If the vectors B, Pi and P 2 are indeed in a 
single plane, then the velocity vectors being proportional 
to BxPi, BxP 2 and Pi xP 2 are all collinear. Let #1 and 
$ 2 be the angles between Pi and P 2 and B, respectively. 
If Pi and P 2 are in the same plane, then both evolution 
equations can be reduced to the form P^ = w c B x Pi 
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(i = 1, 2) with 



sin(t? 2 - gi) 



/ CiS 2 

= +W +|UC( C 2 

sin($ 2 — #i) 



L0 C = -LUQ + fj, ■ 



-UJ a + \1 a 



sin $2 
c 2 si 



S2 



(58) 



where we have used the notation ci = cos $i and so forth. 
This is equivalent to Eq. (65) of Duan et al. [12]. 

Eliminating /j, sin(^ 2 — $i) from the first form of these 
equations we find 



si + as 2 



(59) 



In the alignment case where d\ — d 2 we immediately 
reproduce 



^synch _ 1 + a 

u>o 1 — a 



(60) 



for the synchronization frequency. 

It is also trivial to eliminate lu c from the two equations 
in Eq. (58) and one finds 



_ wo _ Si — as2 I c\ c 2 
~ (j, 2 \si s 2 



(61) 



In the alignment case, #i = i? 2 , the pure precession mode 
requires k — (or \i — > oo). 



They have angles relative to the B-direction of 

Ll)q — U) c + flCl — liOLCl 



cos 6*i 
cos 6*2 



sj (uj - uj c + na - iiaci) 2 + /z 2 (sf - as 2 ) 2 ' 
1/ (-w -w c + /xci - iiac 2 ) 2 + M 2 ( s i - as2) 2 



(64) 



If we insert on the r.h.s. of each of these equations the 
expression for oj c from the corresponding lines in Eq. (58) 
we find that indeed cos #1,2 = cos-^1,2. Therefore, every- 
thing is self-consistent: The pure precession mode is a 
form of motion of both vectors where each is collincar 
with its Hamiltonian. 



C. Lepton-number conservation 

With some mild restrictions as noted above, the pure 
precession mode is possible for almost any combination 
of ■&! and i?2- In view of our application it is useful, how- 
ever, to group this two-dimensional set of solutions into 
families of solutions that have a fixed "lepton number," 
meaning solutions for which 



D\\ = cos^i 



! cos ^2 



(65) 



is a common fixed number. Each of these families in- 
cludes the alignment case $1 = ^2 = ^0 where D\\ = 
(1 — a) cos??o- Therefore, we can write the condition of 
lepton-number conservation in the form 



(1 — a) cos??o = cosi?i — a cos $2 • 



(66) 



B. Alignment with Hamiltonians 

It is surprising that for any two angles $1 and $2 one 
finds solutions lu c and [i that permit a pure precession. 
Of course, suitable initial conditions are required to ac- 
tually put the system into this mode. Demanding [i to be 
positive (as it would for neutrino-neutrino interactions) 
implies some restrictions on the allowed angular ranges, 
but otherwise all combinations of angles are possible. 

As a check of self-consistency we note that in the pure 
precession mode, Pi^ 2 do not move within the co-rotating 
frame. Therefore, in this frame their Hamiltonians must 
be collinear with Pi, 2- We can write the EOMs in general 
as 

Pi = Hi x Pi, 

P 2 = H 2 x P 2 . (62) 

The Hamiltonians in the co-rotating frame are 

Hi = (+w - w c )B + ft (Pi - P 2 ) , 

H 2 = (-w - ^c)B + n (Pi - P 2 ) . (63) 



In other words, once a has been chosen, we can classify 
the solutions by the parameter $0 or equivalently by D|| . 
All solutions with the same i?o form a family with the 
same lepton number and thus are adiabatically connected 
to the aligned case where $1 = 1?2 = $o- 



D. Explicit solution 

We are now in a position to find explicit solutions of the 
EOMs in the adiabatic limit where the two polarization 
vectors move in the pure precession mode. To this end 
we specify a value for a and $0, an angle that has the 
interpretation of twice the effective mixing angle. We can 
then use one of the tilt angles $1 or §2 as a parameter to 
characterize different solutions. 

We know that the shorter anti-neutrino polarization 
vector P 2 for /i — > will be anti-aligned with B whereas 
P2 will remain at some finite angle fixed by lepton- 
number conservation. Therefore, $2 is a more useful in- 
dependent parameter. More specifically, we use 

77 = cos $2 (67) 
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as an independent variable that varies between cosf^o > 
7] > — 1. After eliminating i?i with the help of Eq. (66) 
it is trivial to find from Eq. (59) explicitly 



\Jl - (co - ac + ar]) 2 + ay/T 



(68) 



where cq 
find 



w o y/\ — (c — ac + ar?) 2 — a-\/l — rf 

cost? = cos2#£g. Likewise, from Eq. (61) we 



\J\ - (cq — ac + ar;) 2 — a-\/F 



2y/l -r) 2 y/l- (co - aco + cm?) 2 
(cq - ac + a??) \/l - f? 2 



- r\\J 1 - (c - ac + a??) 2 • (69) 

While we cannot extract an explicit analytic result for 
uj c (k), we have an analytic form of this curve as u) c {v) 
and k(t]) in parametric form. Therefore, in contrast to 
Rcf. [12] and in contrast to our box-spectrum example 
we do not rely on a numerical solution: All properties of 
the solution can be understood analytically. 

In Fig. 11 we show u> c (k) as a thin line for Co = vO/75 = 
0.866, corresponding to sin 26^ — 0.5. This is equivalent 
to the thin line in the upper panel of Fig. 9, except that 
here we have two polarization vectors whereas there we 
had a double-box spectrum. The solutions are very sim- 
ilar. Using a smaller and smaller sin 28, corresponding 
to c approaching 1, the solution develops a sharper and 
sharper cusp as expected. We show the limiting curve 
for cq = 1 or sin 2#fS = as a thick solid line. 
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FIG. 11: Common precession frequency for our system of 
two polarization vectors with a = 0.7, in full analogy to the 
double-box spectrum case of Fig. 9 (upper panel). The thick 
line is the limiting solution for sin 26^}- — » 0, equivalent to 
cosi?o — » 1, t he th in line is for sin 20£§ = 0.5, corresponding 
to costfo = 70.75 = 0.866 as in Fig. 9. 



E. Cusp 

The limiting case cos i?o — > 1 involves a number of 
subtleties. In the strictly adiabatic limit this solution 



corresponds to both polarization vectors being strictly 
aligned with B or, in the neutrino language, to an ex- 
actly vanishing effective mixing angle. Starting with this 
initial condition, reducing the neutrino-neutrino interac- 
tion strength \i could never lead to a deviation of the 
polarization vectors from this orientation just as in the 
absence of mixing there can be no flavor conversions. 

We first study this limiting case by using c = 1 in 
Eqs. (68) and (69). The functions u> c (v) and then 
give us only the large-K branch above the cusp. Specifi- 
cally for t) — > 1 we find the familiar cusp values [11, 12] 



(l-v^) 2 



2 

0.01334 

i + 



for a = 0.7, 



Wo 



l-y/a 
= 11.2444 for a = 0.7. 



(70) 



Numerically these results are very similar to the double- 
box spectrum case, but not identical. 

Let us begin with two polarization vectors in the pure 
precession mode with intermediate \i and let us assume 
the conserved D» is such that in the alignment case they 
have c = 1. If we now increase \i adiabatically, align- 
ment with each other and with B is reached for fi = [i cusp 
so that i?i = $2 = 0. As we increase \x further until in- 
finity, they cannot get more aligned. This is different for 
all cases with c < 1 where the angles keep monotoni- 
cally decreasing with increasing /i and reach alignment 
i9i = = only for infinite \x. 



F. Sleeping-top regime 

The range /x cusp < /i < oo in the limiting case has 
the interpretation of the "sleeping top regime" [12]. In 
the language of the gyroscopic flavor pendulum it corre- 
sponds to a spinning top that spins so fast that energy 
and angular momentum conservation prevent it from de- 
viating from perfect vertical alignment if it was pre- 
pared in this state, seemingly defying the force of gravity. 
Therefore, we call this the sleeping-top range of /x or k. 

We can study this regime using a somewhat different 
paramctrization that allows us to take the limit of a small 
angle #o- To this end we write $2 = £$0 an d thus express 
the tilt angles as multiples of d . In Eqs. (68) and (69) we 
thus substitute Co = cos$o and rj = cos $2 = cos(£i?o)- 
This is only a re-paramctrization: Once more we can 
find, for any chosen value of i?o> the solution uj c (k) in the 
parametric form w c (£) and «(£) with 1 < £ < 7r/?V 

In the limit do — > we find for the co-moving frequency 



y/1 + a(e - 1) + a£ 
VI + a(e - 1) - a£ 



(71) 



This leading term in a t9 expansion does not depend on 
$0 and thus applies even for 1)0 = 0. While in this case 
^2 = C^o = 0, this result nevertheless depends on £. 
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For £ = 1 we reproduce the expected synchronization 
value of Eq. (60) whereas for £ — > oo wc find lo cusp of 
Eq. (70). In other words, in the limit i?o —> our new 
parametrization covers the sleeping-top branch, but not 
the region beyond the cusp. 

One can take a similar representation for and use 
only the leading term. This gives the limiting values 
k(0) = and k(oo) = K cusp as expected. 

Actually, Eq. (71) can be inverted to provide £(w c ) and 
the result can be inserted in the leading term of k(£), 
providing after some simplifications 



1 



a 



w c /w + 1 w c /w - 1 ' 



(72) 



This function is well defined in the entire range w C us P 



< 



lo c < oo, i.e., once more we find both a physical and an 
unphysical branch of the solution. We can also invert 
this result and find for the physical branch 



w c _ 1 - a - y/a 2 - 2a(l + 2k) + (1 - 2k) 2 
ojq 2k 



(73) 



whereas the unphysical branch has a positive sign of 
the square root. This function has real values only for 
< k < z^cusp and indeed approaches the familiar lim- 
its a> sync h and uj cusp at the two ends of the sleeping-top 
interval. 



G. Sum rules 

The case of two polarization vectors studied in this 
section corresponds to the spectrum 



P w = 5(oj - lo ) + aS(uj + uq) . 



(74) 



The adiabatic solutions must obey the sum rules dis- 
cussed earlier. For this two-line spectrum Eqs. (30) 
and (31) give 



1 = 



1 



K — TK -\- Du) 2 



a 



{-K-rK + D«) 2 + D\ 



^(K-rK + D {l ) 2 +D 



+ 



a 



^{-K-rK + D^) 2 + Dl 



(75) 



where r = uj c /^o- For any two angles #i and 1? 2 we can 
calculate the corresponding Dm, D±, lu c and k and show 
that the sum rules are actually fulfilled. 

In the general case we used the sum rules to find so- 
lutions lu c (k) and Dj_(«) and then find solutions for the 
individual polarization vectors which here would mean 



to find i?i(k) and ^(k). In the two-polarization vector 
case we were able to proceed in the opposite way and 
find, for a given i? (or a given D||) the solutions #i, oj c 
and function of #2- This was possible because 

for two polarization vectors lepton-number conservation 
uniquely fixes the angle $i as a function of 1? 2 - 



H. Spectral split? 

For two polarization vectors a spectral split cannot oc- 
cur because it would amount to the polarization vector 
Pi breaking apart. The true final state consists of P 2 
indeed anti- aligning with B and Pi retaining a non-zero 
transverse component. The Pi final orientation is de- 
termined by lepton-number conservation. In the end it 
precesses freely around B with frequency ujq. Indeed, we 
have found that in the end u>® = uio- 

At the same time we have noted that in the pure 
precession mode the polarization vectors are aligned 
with their Hamiltonians which in the end seem to be 
(anti)parallel to B because fj, = 0. Notice, however, that 
Hi in Eq. (63) involves (wo — uj c )B and if in the end 
lu c — > luq the final length of Hi approaches zero. In a lim- 
iting sense it always maintains a direction tilted against 
B such that Pi always remains aligned with Hi. 

Therefore, we here have a case where the perfectly adi- 
abatic evolution leads to a final configuration where one 
of the polarization vectors is not (anti)aligncd with B. 

For an arbitrary number of discrete polarization vec- 
tors an exact spectral split will usually not be possible 
because typically one of vectors would have to be bro- 
ken apart to achieve a split. In this case the final ui c 
will be the u of the left-over polarization vector that can 
neither align nor anti-align with B if lepton number is 
to be conserved. Therefore, we expect that in this sit- 
uation the final adiabatic solution will be such that all 
polarization vectors (anti) align with B except for this re- 
maining one that sits on the fence and for which in the 
end (u> — uj c ) — > 0. In this way it remains aligned with 
its Hamiltonian alright. We have seen such a case in the 
fully adiabatic example in Fig. 8 where one mode ended 
up at P u u = 0. 



VII. ADIABATICITY VIOLATION 

A. Sharpness of the spectral split 

In practice, fj,(t) decreases with a finite speed. On the 
other hand, the P w with frequencies close to iv sp ut precess 
very slowly at late times when fj, — ► 0, so that adiabatic- 
ity is broken for those modes as discussed earlier. As 
a consequence, the spectral split is not infinitely sharp. 
There is a transition region of non-zero width. A quan- 
titative estimate of the adiabaticity condition is given in 
our earlier paper [14]. 
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The width of the transition region, i.e., the sharpness 
of the spectral split, is a common measure of the degree 
of adiabaticity also for the ordinary MSW effect where 
the split occurs at u = 0. 

In order to quantify the meaning of "sharpness" we 
consider the example of neutrinos plus antineutrinos with 
a double-box spectrum with the average frequency lo as 
in Sec. V and study the final neutrino spectrum. From 
Fig. 7 it is evident that the (anti) alignment of all po- 
larization vectors with B is almost perfect, even if adi- 
abaticity is violated, except for the modes around the 
split. In Fig. 12 we show the final P\\(lj) and P±(ui) as 
a function of (w — u) sp ut) / Aw . Here, Auj is the variance 
or root-mean square width of the distribution P±(cu). Of 
course, since P^ is conserved both distributions contain 
the same information, but P± lends itself more directly 
to a straightforward definition of the width of the non- 
adiabatic range. 



find that 



-0.5 




-4-2 2 4 

("-"split)/^ 

FIG. 12: Final spectra of P\\(uj) (upper panel) and P±(u>) 
(lower panel) for a double-box spectrum of neutrinos plus an- 
tineutrinos as in Sec. V. Here, Auj is the rms width of the 
distribution which depends on the time scale r for the expo- 
nential fj,(t) variation. 



In this form, the shown functions appear to be uni- 
versal for the assumed double-box spectrum as far as we 
can tell from a numerical exploration of parameters. In 
particular, these functions do not seem to depend on the 
assumed asymmetry parameter a, not on the value of 
t^spiit an d not on the mixing angle. They do not de- 
pend on the degree of adiabaticity, i.e., on the time scale 
t of the exponential fi(t) variation, except through Au>. 
Varying loqt over two decades between 10 and 10 3 we 



Aw w wo 



(w r)3/4 



(76) 



is an excellent representation with an accuracy better 
than a few percent. (For uj t = 100 this implies Auj = 
0.032 ljq.) However, this relation is not exact: The de- 
pendence of Aw on r does not seem to be an exact power 
law. We give some explanations of these features in the 
following subsection. 



B. MSW effect vs. collective transformations 

To shed more light on the dynamics of the split phe- 
nomenon we compare here the adiabaticity violation with 
its counterpart in the usual MSW case where the evolu- 
tion is given by Eq. (4) with fi = and where the EOMs 
for different frequencies uu decouple. 

Let us consider first the MSW case for the normal mass 
hierarchy. If the initial density is very large, much larger 
than the resonance density, A 3> A res (u;), the initial ef- 
fective mixing is strongly suppressed. The adiabatic evo- 
lution to small densities, A — > 0, then leads to the well- 
known strong conversion with a final survival probabil- 
ity P = sin 2 9 s=s 0. The neutrino polarization vector 
initially coincides with the Hamiltonian. Following the 
Hamiltonian completely flips the flavor. With increasing 
energy E (decreasing frequency w) the degree of adia- 
baticity decreases and above a certain energy E na the 
effects of adiabaticity violation increase exponentially. P 
increases and approaches 1. In the antineutrino channel 
for E — > oo (uj — > 0) the mixing is strongly suppressed 
and the survival probability So, there is a continu- 

ity in the w-variable around u> ~ 0, i.e., at the transition 
from neutrinos to antineutrinos. 

The split frequency can be estimated using the 
Landau-Zener probability which describes well the sur- 
vival probability for small mixing angles [22, 23]: P « 
P LZ = e~ Auj . Here A = ttt sin 2 26/(2 cos 2(9), and 



t = A/A is the scale height of the density change. From 
the condition Plz = 1/2 we find 



t^split 



2 In 2 cos 29 
ttt sin 2 29 



(77) 



For the width of the split region one finds Auj ~ w sp nt and 
Aw/u> sp iit is a universal quantity that does not depend on 
the parameters of the problem. According to Eq. (77) the 
sharpness of the split is 



Auj cx — 
r 



(78) 



which differs from the dependence in Eq. (76) for the 
collective transformation case. 

We emphasize that here as well as in the collective 
transformation case the split phenomenon requires a spe- 
cific initial condition: a small effective mixing angle. 
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For the MSW case in the normal mass hierarchy the 
split is in the neutrino channel and the flavor flips for 
u> > w S piit. For the inverted hierarchy the split is in the 
antineutrino channel and the modes with ui < uj sp ut flip. 

With an increasing degree of adiabaticity (increasing 
t) the split frequency shifts to smaller values, approach- 
ing w sp iit = which is realized in the ideally adiabatic 
case. Conversely, with a decreasing degree of adiabaticity 
both w S pnt and the width of the transition region increase. 
This differs from the collective transformation case where 
w S piit does not shift and is determined solely by lepton- 
numbcr conservation. In other words, the changes of u c 
and D± adjust themselves self-consistently in such a way 
that w S piit = const. In the MSW case flavor lepton num- 
ber is not conserved. 

Let us consider the origin of the differences between 
these two cases. The Hamiltonian H u = H w • er, where 
er is the vector of Pauli matrices, has the diagonal and 
off-diagonal elements ±H and H , respectively, 

2H = nD\\ - {u-u c ), 2H = fiD ± . (79) 

The resonance condition, H = 0, can be written as 



LU = ^ ros -D|| + CJ c (^ ros ). 



(80) 



The general form of the adiabaticity condition is c ff <C 
AH, where # e ff is the effective mixing angle and AH the 
difference of eigenvalues. In the resonance region this 
condition takes on the form 



H 



< 1. 



AH 2 

Using Eq. (79) we find explicitly 



< 1. 



r 

The transition to the usual MSW case is given by 



(81) 



(82) 



(u - UJ C ) 
D± 



A, 

L|| = L cos 29 , 
Lj_ = i sin 26*. 



(83) 



A substantial difference is that in the MSW case both 
L|| and L± are constants and determined by the vacuum 
mixing angle. Moreover, oj c is absent. The resonance 
condition in the mass basis is: u> = A ros cos20. The adi- 
abaticity condition in the resonance becomes 



2t x lo sin 2 29 cos 26 



< 1. 



(84) 



An essential difference of the collective case is that D± 
is a dynamical variable and as such depends on fi. There- 
fore the effective mixing angle changes not only due to the 
shrinking of the vector ^.D in analogy to AL in the MSW 



case, but also due to the rotation of D. Consequently the 
width of the resonance layer, where 29 c f{ changes from 
45° to 135°, and the adiabaticity condition essentially do 
not depend on the initial mixing angle. This explains 
the universality of the function that describes the spec- 
tral split. 

The l.h.s. of the inequalities (82) and (84) are adi- 
abatic parameters which, up to numerical coefficients, 
give the parameters in the exponents of the correspond- 
ing Landau-Zener probabilities. They determine how the 
sharpness of the spectral split depends on various param- 
eters. Apparently Eq. (84) reproduces the result Eq. (78). 
In the collective case the dependence of Aw on r M is more 
complicated due to the presence of uo c and D± which de- 
pend on ii. These functions are given in Fig. 9. Ap- 
proximating these functions with power law functions in 
the resonance region and inserting them into Eq. (82) 
can provide an explicit expression for the shape of the 
washed-out spectral split. 



C. Initial alignment of polarization vectors 

When we identified a sequence of static solutions with 
fixed lepton number and variable fj, as the adiabatic solu- 
tion for a slow fj,(t) evolution, we assumed that initially 
fi — > oo. We assumed that initially all neutrinos are pre- 
pared in weak interaction eigenstates, corresponding to 
all P w being aligned in a direction that is tilted with the 
angle 29^ relative to B. This initial condition of perfect 
alignment corresponds to a pure precession solution only 
if /i — > oo. 

In a realistic system /i cannot be infinite. Therefore, 
if we begin with an initial condition of perfect alignment 
combined with a finite ji, the system is initially not in 
a pure precession mode. Even if the subsequent /i(t) 
evolution is infinitely slow, the polarization vectors will 
not be perfectly aligned with the Hamiltonians of the 
adiabatic solution. 

In a linear system this situation would not spell a de- 
viation from adiabaticity. Usually adiabaticity is taken 
to mean that the polarization vectors follow the Hamilto- 
nians, but they need not be aligned with them. Rather, 
the precession cones would retain a non-zero but fixed 
opening angle. 

Here the situation is different in that the P w define 
D and thus affect the single- mode Hamiltonians H w . 
Therefore, these Hamiltonians must now themselves show 
fast motions when they are viewed from the co-rotating 
frame. Still, if the initial \i is large (but finite), the de- 
viation from a static solution is small and one expects 
that nonlinear effects appear only at second order in the 
small deviations of the polarization vectors from their 
adiabatic solutions. One would expect that the P w pro- 
cess around the nearly static H w where the opening an- 
gle of the precession cone no longer vanishes. Notice that 
since D depends on many modes the fast motion of the 
latter to some extend cancels in D and therefore H w . 
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This behavior is borne out in numerical examples. 

This numerical observation implies that the adiabatic 
solution is a stable fixed point of the system, i.e., a small 
perturbation does not lead to a run-away. Rather, the 
system always evolves close to the adiabatic solution and 
oscillates around it. The stability of the adiabatic so- 
lution is not automatically implied by its existence as a 
self-consistent solution of the EOMs. Since we do not 
have developed a formal criterion of stability, we do not 
know if there could be examples that would not be sta- 
ble. Therefore, it is a numerical observation, not a proven 
fact, that the adiabatic solution is a good approximation 
to the exact solution for typical cases of realistic initial 
conditions. 



D. Nutations 

A system consisting of neutrinos and antineutrinos is 
similar to the neutrino-only case if one interprets the 
antineutrinos as negative-frequency modes, but it has a 
number of peculiar properties. One such a feature con- 
sists of the "wiggles" in the P u ii evolution that we noted 
in the numerical solution of the EOMs shown in Fig. 8. 
The amplitude of these modulations becomes deeper if 
the effective mixing angle is smaller. To illustrate this 
point we show in Fig. 13 (top panel) the same example 
as in the upper-left panel of Fig. 8, now with a smaller 
mixing angle sin2#^j = 1CP 3 . 

Our general treatment shows that in the limit of a very 
slow ^(t) variation and if initially fi — > oo all polarization 
vectors move in a co-rotating plane where the EOMs are 
solved by our adiabatic solution without wiggles. They 
are, therefore, a non-adiabatic feature and indeed disap- 
pear numerically if we use an unrealistically slow fi(t) 
variation. 

It has been shown that the EOMs of a system con- 
sisting of only two polarization vectors, one for neutrinos 
and one for anti-neutrinos as in Sec. VI, are equivalent 
to a gyroscopic pendulum in flavor space [11, 16]. Its 
general motion consists of a precession around B and 
a nutation that causes the overall zenith-angle modula- 
tion. For very large neutrino-neutrino interactions (very 
small k), the pendulum is in the "sleeping top phase" 
[16] where its motion is a nearly pure precession. In the 
neutrino language this corresponds to the synchronized 
regime where all polarization vectors are strictly pinned 
to each other and precess around B as a collective object. 
Once k becomes larger than the "cusp value," the gyro- 
scope begins to wobble. In the neutrino language, this is 
the point where the polarization vectors begin to spread 
in the zenith-angle direction and arc no longer strictly 
pinned to each other. 

The nutations are excited because the flavor pendulum 
initially deviates only by the small angle 29^ from the 
inverted vertical direction. Starting the pendular motion 
requires this small angle to grow under the influence of 
the external force, introducing a time scale that depends 




FIG. 13: Top: P^||(/i) for 25 individual modes for the case 
of neutrinos plus antineutrinos. Same as upper left panel of 
Fig. 8, but now with a small mixing angle sin 26^ = 10~ 3 . 
Middle: D± calculated from the numerical solution compared 
with the adiabatic solution (thick line). Bottom: M±/D± in 
units of u>o calculated from the numerical solution compared 
with the adiabatic lo c (thick line). 



logarithmically on 29^. If the fi(t) variation is fast com- 
pared to this time scale, the pendulum cannot follow the 
adiabatic solution. This explains that the depth of the 
nutation amplitude becomes larger for a smaller 29^. 
The case of a vanishing 29^ was the main example that 
we studied in Sec. V. While one can find a perfectly self- 
consistent static solution for this case, it is clear that the 
system cannot follow this solution: For a vanishing mix- 
ing angle the system will not move at all, i.e., the delay 
leading to nutations is infinite. Therefore, this example 
was an instructive limiting case, but cannot be realized as 
an adiabatically connected sequence of static solutions. 

A remarkable feature of the nutations shown in Fig. 13 
is that they are a collective effect of all modes. We cannot 
interpret the nutations as a simple precession of the indi- 
vidual P w around the static Hamiltonians. Rather, the 
nutations represent a motion of the entire system around 
the adiabatic solution. To illustrate this point we show 
in the middle panel of Fig. 13 the evolution of the numer- 
ical D± in comparison with the corresponding adiabatic 
result. The true solution oscillates around the adiabatic 
solution with the nutation period. Likewise, the instanta- 
neous precession frequency of the co-rotating plane, given 
by the quantity M±/D±, oscillates around the adiabatic 
result for lu c (bottom panel). 

Once more we are numerically led to the conclusion 
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that the adiabatic solution represents a stable attractor 
of the nonlinear system. The true solution is a collec- 
tive oscillation mode around this solution and thus plays 
the role of a perturbation. At the same time we stress 
that the individual P w evolve under the influence of their 
instantaneous Hamiltonians H w which are not identical 
with the ones of the adiabatic solution. The instanta- 
neous Hamiltonians themselves nutate, although they nu- 
tate around the static Hamiltonians. We recall, however, 
that the instantaneous Hamiltonians always lie in a sin- 
gle plane, whereas this is not the case for the polariza- 
tion vectors. Moreover, the motion of the instantaneous 
Hamiltonians is governed by the variation of the single 
vector D so that they must move in lockstep with each 
other. 

Another important feature of the nutations is that they 
are a transient phenomenon that dies out toward the end 
of the evolution when /i — > as explained in the flavor 
pendulum picture [11, 16]. As a consequence, the spectral 
split remains unaffected. Even though all polarization 
vectors wildly nutate, they begin and end (anti) aligned 
with B except for those modes near the split that anyway 
do not reach the final (anti) alignment. The modes suffi- 
ciently far from the split begin and end aligned with the 
static Hamiltonians, whereas for intermediate times they 
wildly nutate. The nutations have no apparent impact 
on the final flavor spectrum of the system. 



VIII. CONCLUSION 

We have elaborated our previous theory [14] of spec- 
tral splits that occur as a consequence of the adiabatic 
evolution of a neutrino ensemble. We have clarified that 
the spectral-split phenomenon in itself is not an entirely 
new effect, but a very useful way of looking even at the 
usual MSW effect if we express its energy dependence 
in terms of the oscillation frequency u) — Am 2 /2E and 
if we identify antineutrino modes with negative frequen- 
cies. However, in the case of collective transformations 
the effect has some distinctive features: the position of 
the split is fixed by the flavor lepton number of the sys- 
tem and docs not depend on other parameters. 

The main impact of collective neutrino-neutrino inter- 
actions is that the evolution is adiabatic not in the lab- 
oratory frame, but rather in a co-rotating frame that is 
defined by a plane in flavor space which is spanned by 
the mass direction B and the vector D which is the po- 
larization vector of the overall flavor lepton number of 
the system. This frame defines a static solution of the 
EOMs that consists, in the laboratory frame, of a pure 
precession of all polarization vectors around B. If the 
evolution in the co-rotating frame is adiabatic, we will 
have a sharp spectral split at the frequency lo — in the 
co-rotating frame, corresponding to a frequency w S piit in 
the laboratory frame where usually w sp ii t ^ 0. 

In order to take advantage of this picture we must find 
the co-rotating frame explicitly This is achieved in terms 



of two sum rules (w-integrals over the spectrum P w ) that 
guarantee the self-consistency of the static solution, i.e., 
of the pure precession mode. 

The main part of our paper is devoted to explicit an- 
alytical and numerical solutions of these equations (sum 
rules) for generic examples. The adiabatic solution is de- 
termined completely by two functions w c (^) and D±(jj,). 
We study properties of these functions for several differ- 
ent neutrino spectra. For the neutrino-only box spec- 
trum and non-zero initial effective mixing angle both 
functions decrease monotonically with [i: uj c (^l) — > a> sp iit 
and L> ± (^) ^ 0. 

In the presence of antineutrinos the behavior substan- 
tially changes: For small initial mixing angles the depen- 
dence Lv c (p) has a cusp at a certain density ^i C usp- 

The 

cusp corresponds to the border between the sleeping-top 
region and the beginning of bipolar transformations. Es- 
sentially the flavor evolution starts for densities below the 
cusp. Above the cusp, w c (/x) increases starting from the 
synchronization frequency. Below the cusp it decreases 
monotonically. With decreasing [i we find that D±(/j,) 
first increases, reaches a maximum at [i <~ uj , a typical 
frequency of the spectrum, and then decreases. 

This type of behavior appears to be universal and 
shows up for completely different spectra. Very similar 
quantitative features appear for the continuous double- 
box spectrum and for the bi-chromatic spectrum. The 
latter case, by the gyroscopic flavor-pendulum analogy, 
allows one to understand a number of subtle features 
which also show up for continuous spectra. 

The adiabatic solution is an approximate solution of 
the EOMs. There are several features which appear in 
the exact solution that are not described by the adiabatic 
solution. The violation of adiabaticity always takes place 
at late stages of evolution (at small /i) for modes with u> 
close to w S piit. This leads to a partial wash-out of the 
split so that the degree of adiabaticity determines the 
sharpness of the split. Numerically it appears that for 
a given spectrum the profile of the split is a universal 
function of the density gradient and does not seem to 
depend on other parameters. Measuring the width of the 
split region would allow one to determine the neutrino 
density gradient in the region /i ~ ljq . 

Nutations related to bipolar transitions represent an- 
other deviation from adiabaticity. The nutations occur 
around the adiabatic solution, or, in other words, the 
adiabatic solution appears as an average over nutations. 
Somewhat surprisingly, the nutations do not visibly influ- 
ence the final result. They disappear as the neutrino den- 
sity decreases and thus are only a transient phenomenon. 

Our treatment leaves a number of interesting questions 
unresolved. While the self-consistency conditions deter- 
mine the static solution, we have not realized under which 
circumstances such solutions actually exist, although we 
have not found a counter example. 

The adiabatic solution is stable in typical numerical ex- 
amples in that the true solution of the EOMs tracks this 
adiabatic solution. There could be situations where this 



20 



is not the case, leading to kinematical decoherence rather 
than a collective oscillation of the individual modes. 

Even though several important aspects of this sys- 
tem are only numerically observed and not analytically 
proven, we believe that our adiabatic treatment of the 
EOMs nicely explains many of the intriguing features of 
this nonlinear system. 
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